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Abstract 

The Chern-Simons bosonization with U{1) x SU{2) gauge field is apphed to 
the 2-D t — J model in the limit t S> J, to study the normal state properties 
of underdoped cuprate superconductors. We prove the existence of an upper 
bound on the partition function for holons in a spinon background, and we 
find the optimal spinon configuration saturating the upper bound on average 
- a coexisting flux phase and s + id-like RVB state. After neglecting the 
feedback of holon fluctuations on the U{1) held B and spinon fluctuations on 
the SU{2) field V , the holon field is a fermion and the spinon field is a hard- 
core boson. Within this approximation we show that the B field produces a 
TT flux phase for the holons, converting them into Dirac-like fermions, while 
the V field, taking into account the feedback of holons produces a gap for 
the spinous vanishing in the zero doping limit. The nonlinear cr-model with 
a mass term describes the crossover from the short-ranged antiferromagnetic 
(AF) state in doped samples to long range AF order in reference compounds. 
Moreover, we derive a low-energy effective action in terms of spinous, holons 
and a self-generated U{1) gauge field. Neglecting the gauge fluctuations, the 
holons are described by the Fermi liquid theory with a Fermi surface consisting 
of 4 "half-pockets" centered at (±^,ib|) and one reproduces the results for 
the electron spectral function obtained in the mean fleld approximation, in 
agreement with the photoemission data on underdoped cuprates. The gauge 
fluctuations are not confining due to coupling to holons, but nevertheless 
yield an attractive interaction between spinous and holons leading to a bound 
state with electron quantum numbers. The renormalisation effects due to 
gauge fluctuations give rise to non-Fermi liquid behaviour for the composite 
electron, in certain temperature range showing the linear in T resistivity. This 
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formalism provides a new interpretation of the spin gap in the underdoped 
superconductors (mainly due to the short-ranged AF order) and predicts that 
the minimal gap for the physical electron is proportional to the square root of 
the doping concentration. Therefore the gap does not vanish in any direction. 
All these predictions can be checked explicitly in experiment. 

PACS Numbers: 71.10.Pm, 71.27.+a, 11.15.-q 
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Typeset using REVT^ 



I. INTRODUCTION 



A. Physical issue to be addressed 

The proximity of superconductivity (SC) to antiferromagnetism (AF) in reference com- 
pounds is a distinct feature of the high-T^ superconductors. Upon doping the AF goes 
away, giving rise to SC. At the same time, the Fermi surface (FS) is beheved to develop 
from small pockets around (±f , ±f ) anticipated for a doped Mott insulator, to a large 
one around (7r,7r), expected from the electronic structure calculations p[ and confirmed by 
the angle-resolved photoemission spectroscopy (ARPES) experiments 0. To understand 
this crossover is one of the key issues in resolving the high Tc puzzle. For this reason, the 
underdoped samples present particular interest due to the strong interplay of SC with AF. 

There is a consensus now that these systems are strongly anisotropic, and the fundamen- 
tal issue is to understand the behavior of strongly correlated electrons in the copper-oxygen 
plane A "spin gap" or "pseudogap" has been invoked to explain the reduction of 
magnetic susceptibility x below certain characteristic temperature T* and suppression 
of the specific heat compared with the linear T behavior This gap also shows up in 
transport properties [|], neutron scattering 0, and NMR relaxation rate measurements. 
The recent ARPES experiments [|lOl on underdoped samples seem to indicate that the FS 
of these compounds is probably half-pocket-like, i.e. a small pocket near S-point (|, |), but 
lacking its outer part in the reduced Brillioun zone scheme. These data show clear Fermi 
level crossing in the (0, 0) to (tt, tc) direction, but no such crossing was detected in the 
(0, vr) to (vr, tt) direction JlOl. The observed pseudogap above Tc is consistent with d-wave 
symmetry. 

In this paper we will be concerned with the normal state properties of these underdoped 
cuprate superconductors, focusing on the implications derived from the proximity of these 
systems to the AF reference state. 



B. A brief survey of related theoretical approaches 

Theoretically there have been mainly two competing approaches: One starting from the 
Mott-Hubbard insulator, advocated by Anderson using the concept of spin liquid, or 



the Resonant Valence Bond (RVB) state, while the other starting from the more conventional 
Fermi liquid (FL) point of view. 

One of the crucial concepts within the first approach is "spin-charge separation" which 
can be intuitively implemented by introducing "slave bosons" ||T^ , namely, one rewrites the 
electron operator: 

where is a charged spinless (slave) boson operator (ho Ion), while fia is a neutral, spin 1/2 
fermion operator (spinon) satisfying contraint 

eUi + flU = l. (1.1) 

(Hereafter the repeated spin indices are summed over.) One can also interchange the role 
of boson and fermion operators, i.e., to introduce a spinless fermion to describe the charge 
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degree of freedom, while "spinning" 



bosons to describe the spin degree of freedom. This 
is the "slave fermion" approach [13,14|. The essential requirement for both approaches is 
the "single occupancy" constraint which is very difficult to implement. In the mean field 
approximation (MFA) which satisfies the constraint ( |1.1D only on average, the "slave boson" 
rSH and "slave fermion" p!^JT6|] approaches gave very different phase diagrams and each of 



them has its own difficulties There have been several attempts to improve the situation 
P^ , but the basic difficulty still remains. 

Moreover, in decomposing the physical electron into a product of fermion and boson, 
one increases the number of degrees of freedom (d.o.f.) by two. The constraint ( [L.ID takes 
care of one, but there is one extra d.o.f. which corresponds to the spinon-holon gauge field. 
In fact, the physical electron operator is invariant under the transformation: 



fjcr — ^ /jo-e 



SO one can "gauge-fix" according to the choice. This is the starting point of the gauge 
field approach to the strongly correlated electron systems |]T^^. This approach has been 
systematically pursued by P.A. Lee and his collaborators first as the U{1) gauge theory ^1 
and recently by considering the SU (2) gauge group 



From the PL point of view the interplay of AF with SC, and the evolvement of the FS 
with doping has been elaborated by Kampt and Schrieffer |2^, and by Chubukov and his 
collaborators |2i|. Very recently, Zhang [^| has proposed an interesting 5*0(5) model to 
consider the AF-SC interplay. 



C. Basic Idea of Chern-Simons Bosonization 



In this paper we will use the C-S bosonization as the basic technical tool. The procedure 
of reformulating the fermion problem in terms of bosons was pioneered by the Jordan- Wigner 



transformation 26 



(1.2) 



where c, is a fermion operator, while a, is a hard-core boson operator on a linear chain. 



The abelian bosonization for one-dimensional fermions with linear dispersion is similar 
in spirit and has found extensive applications in condensed matter physics p8| , p9 | 
formula there is the Mandelstam representation (in field-theoretical jargons) 



The key 
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^^J^^dy^iy) - (-l)"^0F0(x) 



1.3) 



where ip is the fermion operator, while is the boson operator, and a = 1,2. (Here we 
have omitted the normal ordering of the operators.) A rigorous derivation of ( |1.3|) in terms 
of path integrals was given in Ref. [PT|. This bosonization procedure was generalized by 
Witten to the non-abelian case 
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which reformulates the 1+1 dimensional relativistic 
fermion problem with symmetry group G as G-valued nonlinear a (NLo") model. That 
scheme has been extensively used for the study of quantum spin systems [Bl]. 
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The abelian bosonization procedure has been generahzed to 2+1 dimensional systems 
34 1 . It is in some sense analogous to the Jordan- Wigner transformation. The typical relation 



IS given as: 



^Pix) ~ 0(x)e^ ^7. My)dy^^ (14) 



where ip{x) is the fermion field operator, while is the boson field operator, and the two 
are related by the Chern-Simons (C-S) U{1) gauge field operator A^. The integration in the 
exponent is taken over an arbitrary path '-fx in the 2D plane, running from x to oo. The 
path integral will contain an extra factor e"^^"-"- with the C-S action 



Sc.s. — j ^^^^^ ^ ^nd,yAp, 



and the C-S coefficient (level) = 1/(2/ + 1), / = 0, 1, 2.... For non-relativistic fermions the 
C-S constraint can be solved explicitly, and the transformation ( p. .41 ) becomes p^ : 



,.2_„2 



where Q{x — y) = arctan|j— |r. An interesting application of this formula is the analysis of 
the fractional quantum Hall effect at filling z/ = 1/(2/ + 1) in terms of boson liquids p5[ . 



The statistical transmutation, implemented by the abelian C-S gauge field is a consequence 
of the Aharonov-Bohm effect, and it is limited to abelian fractional statistics, characterized 
by the phase factor 2n9 with 9 G [0, 1) In particular, 9 = 0,^ corresponds to boson and 
fermion, respectively, while 9 = j corresponds to the "semion" case, advocated by Laughlin, 
as a constituent quasiparticle in high temperature superconductors ||3^ . 



The 2D analogue of the Jordan- Wigner formula ( |1.4| ) and ( |1.5| ) was originally derived as 
an operator identity [^]. It was later on justified in the path integral form and generalized 



to the non-abelian case [Q. That paper was further extended to include the correlation 



functions, elaborated for the case G = U{1) x SU{2) and applied to the t — J model ||39 



see also PD|). Readers are referred to those references for a detailed presentation. To 



make the present paper more self-contained, we briefiy outline here the basic idea of such a 
bosonization procedure. 

Consider a system of N spin | fermions ( or bosons ) in two space dimensions, in an 
external (abelian) gauge field A. The canonical partition function in the first quantized 
path integral representation is given by 



where aj = 1, 2 for j = 1, 2, ....A^ are the spin indices, a;j(r) represent "virtual trajectories" of 
particles going from imaginary time r = 0tor = /5 = l/T, with T as the temperature (we set 
the Boltzmann constant fc^ = 1) and reaching the plane t = (3 ai the same set of points where 
they start at r = 0, the points being arbitrarily permuted. Due to the fermion statistics of 
the particles, there is a factor (— 1)'^*^''^ associated with each permutation vr, where a{Tc) is the 
number of exchanges in the permutation vr. These trajectories have vanishing probability 
to intersect each other for a given r. Each set of such trajectories appearing in the first 
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quantized representation form a link, i.e. a set of possibly interlaced loops, when the r = 
and T = j3 planes are identified by periodicity in time. Here J^A = J {Aidu^ + Aodr) is a 
line integral in 2+1 dimensions, I = 1,2. The partition function for the boson system would 
be the same, except that a{7i) is replaced with 0. On the other hand, the factor (— 1)'^('^) is 
a topological invariant naturally associated with the link and, according to a general theory 
P^ , it can be represented as the expectation value of a "Wilson loop" (trace of a gauge 
phase factor) supported on that link in a gauge theory with a suitably chosen C-S action. 



where is a C-S gauge field with symmetry group G, k is the C-S cofficient (level), already 
defined in the abelian case, P(-) is the path-ordering, identical to the time-ordering, if "time" 
is parameterizing the path, and 

Scs. = 4^ / d'xTr[e''''f{V,d,V, + '^V.V^V,)]. (1.7) 

As a consequence, there is a boson-fermion relation for the canonical partition function: 

Zp{A) = J DVZB{A + V)e~^^'^-'-^^\ 

The bosonization formula is written in the second-quantized form for the grand-canonical 
partition function: 

^ = E 

and 

E,{A) = Jd^ Z^M/*e-^(^' ^) = / D$ D$* DVe-[^(^' <^*,A + V) + kSUV)] , 

(1.8) 

where \1/ is the Grassmann variable representing the fermion field, while $ is the complex 
variable representing the boson field. 

There might be different choices of G, k: In particular, G = U{1), k = 1 corresponds 



to the abelian C-S bosonization [Q, while G = U{1) x U{1), k = {6, corresponds to 



the "anyon" bosonization ^3[, with 9 as the statistics parameter. In this paper we will 
concentrate on the case G = U{1) x SU{2) with = 2, 1, respectively 
The relations for the correlation functions are given by 

m{x) P{e ^y-^)^{x), ^*(x) ^ $*(x)P(e~*^7. ^) (1.9) 

with 7^. as a straight line in the fixed time plane joining x with oo. This is a non-abelian 
generalization of the 2D formula for the abelian case (|1.4| ). It is important to notice the 
non-local character of these relations. 

Here we have used the first- quantized form of the path integral to identify the relation 
between the fermion and boson systems, while the "working" formula for bosonization is 
given in the second-quantized form of the path integral representation. This switching from 
first to second quantized form and vice versa will be frequently used throughout this paper. 
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D. Outline of the Paper Content 



The f/(l) X SU{2) C-S bosonization approach has been successfully employed by us ear- 
lier to calculate the critical exponents of the correlation functions in the ID t — J model 
in the limit t ^ J. Although, in principle, all bosonization schemes should yield an exact 
identity between the correlation functions of the original fermionic field and correspond- 
ing bosonic correlation functions, the MFA, as mentioned earlier, gives different results in 
different bosonization schemes. The U{1) C-S bosonization has been shown to correspond 
essentially to the slave-boson and slave-fermion approaches (depending on the choice of the 
gauge fixing); while the non-abelian U{1) x SU (2) C-S bosonization corresponds to the slave 
semion-approach {9 = i) ||38|. We have shown |^ that the "semion" spin-charge separation 



of spinon and holon is the correct one to reproduce the exact exponents known from the 
Bethe ansatz solution and the Luttinger liquid-conformal field theory calculations . 



We considered ID fermion system on the background of a 2D U{1) x SU{2) C-S gauge 
field. The U{1) field is related to the charge, while the SU{2) field is related to the spin 
degrees of freedom. Performing the dimensional reduction and using the freedom in gauge- 
fixing, we could analyse the original problem as an optimization process for the partition 
function of holons in a spinon background. We could find an upper bound for the partition 
function and an exact way to saturate this bound, without any approximations. Afterwards 
we used MFA to consider an "averaged" holon configuration to compute the long-time, large 
distance behavior of the correlation functions, reproducing the exact results. The important 
lesson we learned there is that the statistical (semionic) properties of the constituent particles 
(spinous and holons) due to gauge field fluctuations are crucial. In the operator form, the 
original fermion operator can be decomposed as: 

i^x = h,e [e 6J (1.10) 

where is a fermion, while bx is a hard core boson operator. However, only together with 
the attached "holon string" (the exponentiated operator) represents a "physical" charged, 
spinless holon, whereas bx together with the attached "spinon string" corresponds to the 
neutral, spin | spinon. Both of them satisfy the "semionic" equal-time commutation re- 
lataions: 

fix')f{y') = e±^/(yi)/(x^), ^ y\ 

where / represents fermion (boson) operator along with the attached string. This means 
that the "semionic" spinous and holons are "deconfined" , and the spin and charge are fully 
separated. The dynamics of the spinon filed bx is described by an 0(3) NLcr-model. Due to 
the presence of the topological term in ID [2S,2^, the spinous are massless, and they are 



"deconfined", i.e., the spinous themselves, not their bound states (the ususal spin waves in 
higher dimensions) are the constituent quasiparticles. These results encourage us to explore 
the 2D case which is of much more physical importance. We should, however, carefully 
distinguish which are the generic features of the C-S gauge field theory under consideration, 
and what is specific for the ID case. 

In this paper we employ the ?7(1) x SU{2) C-S bosonization scheme to study the 2D 
t — J model in the underdoped regime in the limit t ^ J. We will try to follow as much 
as possible the same procedure as in ID. The U{1) gauge field B is again related to the 
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charge degree of freedom, while the SU{2) gauge field V is related to the spin degree of 
freedom. First we prove the existence of an upper bound of the partition function for holons 
in a spinon background, and we find the optimal, holon- dependent spinon configuration 
which saturates the upper bound in an average sense. The optimization arguments suggest 
coexistence of a flux-phase with an s+zd-like RVB state where the expectation value 
for the AfHeck-Marston(AM) bond-variable of spinous is close to 1, while the s + Z(i-RVB 
order parameter is much smaller than 1. Then we make an approximation, neglecting the 
feedback of holon fluctuations on the U{1) field B and spinon fluctuations on the SU{2) field 
V. Hence the holon field is a fermion and the spinon field a hard-core boson. Within this 
approximation we show that the B field produces a vr-flux phase for the holons, converting 
them into Dirac-like fermions, while the V field, taking into account the feedback of holons 
produces a gap for the spinous, minimal at (±|,±|). The spinous are described by a 
NLcr-model with a mass term (gap) ~ a/— 5 In 5, with 6 as the doping concentration. This 
corresponds to a short-ranged AF order (or disordered state in the jargons of NLcr-model) 
in doped samples, which crosses over to the long-ranged AF order in the prestine samples, 
when the gap vanishes. To our knowledge, this is the first successful attempt to include 
AF fluctuations self-consistently in the RVB-type approach. Moreover, we derive a low- 
energy effective action in terms of spinous, holons and a self-generated U{1) gauge field. 
Neglecting the gauge fluctuations, the holons are described by the FL theory with a FS 
consisting of 4 "half-pockets" centered at (±|,±|) and one reproduces the results for the 
electron spectral function obtained in MFA (for the co-existing 7r-flux and rf-wave RVB 
state) in qualitative agreement with the ARPES data |[T^ for underdoped cuprates. If 
the gauge field were coupled to the spinous alone, it would be confining, since the spinous 
are massive. However, due to coupling to the massless branch of holons (which are actually 
non-relativistic because of a finite FS), gauge fluctuations are not confining, but nevertheless 
yield an attractive interaction between spinous and holons leading to a bound state in 2D 
with electron quantum numbers. This could explain why neglecting the "semionic" nature of 
spinous and holons is less dangerous in 2D than in ID. This means that the spin and charge 
are not fully separated (like in ID), showing up as bound states in low-energy phenomena. 
The renormalisation effects due to gauge fluctuations would lead to non-FL behaviour for 
the composite electron, in certain temperature range showing the linear in T resistivity 
discussed earlier This formalism provides a new interpretation of the spin gap in 

the underdoped superconductors (mainly due to the short-ranged AF order) and predicts 
that the minimal gap for the physical electron is proportional to the square root of the doping 
concentration. Therefore the gap does not vanish in any direction. All these predictions can 
be checked explicitly in experiment. 

The C-S gauge field approach has also been used by Mavromatos and his collaborators 
^ to study the anyon superconductivity, advocated by Laughlin [^]. To our understand- 



ing, the basic aim of their work is to construct a model exhibiting semion superconductivity 
without breaking the time-reversal and parity symmetry. In spite of some apparent similar- 
ities in formulas, the main issue considered and the basic physical assumptions in their work 



are very different from ours. They have also discussed the normal state properties [g^, but 



the mechanism for a possible Non-FL behaviour in their paper differs from what we consider 
here. We should also mention that the SU{2) gauge field considered in their recent papers 
(quoted in pT|) corresponds to a generalization of the local SU{2) symmetry at half- filling, 
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and is not related to the spin rotational symmetry, as we discuss in this paper. P.A. Lee and 
his collaborators |22] have also been considering this (rather than the spin) SU{2) group. 

The present paper is an extended version of the earlier short communication p3[ . 

The rest of the paper is organized as follows: 

In Sec. II we summarize the f/(l) x SU (2) bosonization in the context of 2D t — J model; 
In Sec. Ill we present the optimization problem for the spinon configuration; 
In Sec. IV we derive the spinon effective action; 
In Sec. V we consider the holon effective action; 
In Sec. VI we make some concluding remarks. 

The proof of the bound employed in Sec. Ill, is deferred to the Appendix. 



II. U{1) X SU{2) CHERN SIMONS BOSONIZATION OF THE t - J MODEL 

A. The Model Hamiltonian 

It is widely believed that the 2D t — J model captures the essential physical properties 
of the Cu — O planes characterizing a large class of high-T^ superconductors The 
hamiltonian of the model is given by 

H = Pg[Y. -tii^l^ja + h.c.) + MJ^i^iP ■ i^n'^M Pg, (2.1) 

<ij> 

where aa/3 are the Pauli matrices, ipia is the annihilation operator of a spin | electron on site 
i of a square lattice, corresponding to creating a hole on the Cu site, and Pq is the Gutzwiller 
projection eliminating double occupation, modelling the strong on-site Coulomb repulsion. 
To simplify notations, we introduce a two-component spinor ipi = {ipi^,ipii). Throughout 
this paper, the small letters will denote operators, while the capital letters will denote the 
corresponding complex ( or Grassmann) variables in the path integral representation, unless 
otherwise specified. 

Using the Hubbard- Stratonovich transformation to introduce a complex gauge field 
X<jj>, the grand-canonical partition function of the t — J model at temperature T = 1/ (5 
and chemical potential /i can be rewritten as [^]: 

St_j(/?, /i) = / VXVX*Vmvm*e~^^~J^^ ^ ^' (2.2) 

with 



S,_j{^,^\X,X*)= / dx'lY. X* X<,,> + [(-t + X<,,>)v[/*„vl>^„ + /^.c. 



+ E ^L(50 + + E ^«,.*L*;/3*.73*.a , (2.3) 



where the two-body potential is given by 



9 




1 = J 

Ui,j = { -J hi n.n. (2.4) 
otherwise. 

Hereafter we denote the euchdean-time x° = r ( 9o = d/dr) and its dependence of the 
fields is not exphcitly spelled out. 



B. C-S Bosonization 



Comparing ( U^) , ( P-3D with (|1.8|) , we find that the C-S bosonization procedure, briefly 
introduced in Sec. can be applied to rewrite the grand-canonical partition function (^.21). 
However, there is an important difference, namely, the consideration in Sec. was for 2D 
continuum, where the probability for two world-lines (Brownian paths) to intersect each 
other at a given time is zero. This is not true for the lattice case we consider now, where the 
probability for two paths to cross each other at a given time is not vanishing. On the other 
hand, the model we consider contains a single-occupancy constraint, expressed in terms of 
the Gutzwiller projection operator Pq, or the infinite on-site repulsion Ui^i, which excludes 
the intersection of paths. Therefore, we can still apply the C-S bosonization scheme, leaving 
the C-S gauge fields in the continuum, while considering the matter field on a discrete lattice. 
We will introduce an abelian U{1) gauge field B related to the charge degree of freedom and 
a SU (2) gauge field V related to the spin degrees of freedom. The euclidean C-S actions for 
these fields are given by: 



Sr. 



{V) = ^Jd^xTrU'^PiV,d,V, + IV,KV,] 



(2.5) 



where = a = 1, 2, 3, /i = 0, 1, 2 with cxa as Pauli matrices. 



In the fermion-boson transformation formula (|1.8| ) k = 2 for the U{1) field B, and k = 1 

The correlation functions for the Grassmann fields \I/q, (^'^) 



for the SU{2) field V |8|,|9 



are substituted by the correlation functions of the gauge- invariant complex fields ^ailx) 



($;(7^)), defined as: 



1^ [Be 



$a(7x) 

As mentioned earlier, 7^. is a straight line in the fixed-time plane joining point x with 

and P is the path-ordering operator. In 



00 (reaching a compensating current at 00 p8|j3^ 
principle, we can choose other gauge groups G to implement the bosonization scheme, but 
the encouraging result for the ID t — J model, reproducing the known exact exponents of 



the correlation functions strongly favours the U{1) x SU{2) choice. 

The bosonized action is obtained via substituting the time derivative by the covariant 
time derivative and the spatial lattice derivative by the covariant spatial lattice derivative 
in the U{1) x SU{2) bosonization 



(9o + 25o(j))l + ^K)(j) 



a/3 
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We now decompose the bosonic field as follows: 

^xa = EJl^a, ^*xa = -^x^L) (2-6) 

where 1 is a unity matrix, E a complex scalar field and a complex two-component spin 
i field satisfying the constraint 

Ka^xa = 1. (2.7) 



The gauge ambiguity involved in the decomposition ( p.6|) will be discussed later. The 
field E is coupled to the U{1) gauge field B and it describes the charge degrees of freedom 
and the field Sq, is coupled to the SU{2) gauge field V and it describes the spin degrees 
of freedom of the original fermion. In this description the nature of the groups associated 
with charge (U{1)^ and spin (^SU{2fj is explicitly exhibited and the coefficients of the C-S 
actions are such that the charged and spin ^ field operators reconstructed from the gauge 

invariant (euclidean) fields [ExC^-i^^ and P{e ^-i^'^)a(3^x(3) obey semionic statistics |]38| , ^ . 
In terms of E and E the U{1) x 5'f/(2)-bosonized action of the t — J model is given by 

St-j{E,E\T.,i:\X,X\B,V) 



+ J2uijE*EiE*Ej} + 2S,.s.{B) + 5,.,.(V^) (2. 



with constraint ( |2.7| ) and Coulomb gauge-fixing for the U{1) x SU{2) field implemented 

It is convenient to describe the charge properties in terms of a hole-like field H (holon) 
and this can be achieved in this formalism by substitution E H*,E* — > H, with H,H* 



as Grassmann fields, and changing the sign of the C-S action for the B field p9|j4^ . After 
integration over the auxiliary gauge field X, the grand-canonical partition function 
can be rewritten as: 

IJ') = J VHVH*V^J)T.lVBVVe~^^^^ ^' ^' ^)(5(S*S - 1), (2.9) 
where the euclidean action is given by: 

S{H, H\ S, S*, 5, V) = dx'{Y^ [H*[do - iBoij) - (/x + ^))h, + tBoU) 
+ (1 - H*H,)^,{do + tVoU)) + E [i-tH;e'^<'^>''H,i:UPe'^<'^>'').pJ:,f, + h.c.) 

<ij> 



In what follows we will denote the shifted chemical potential fi' = fi + J/2 hy 6, proportional 
to the doping concentration. 
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C. Gauge Fixings 



The action ( |2.1CI| ) is invariant under the local gauge transformations: 

f/(l) : Hj Hje'^^^\ H* H*E-'^^j\ 

B^{x) ^ B^{x) + d^A{x), A(x) e R 

SU{2): ^^^R\j)T.^, S*^S*i?(j), 

V^,{x) R\x)V^{x)R{x) + R^{x)df,R(x), R{x) e SU{2) (2.11) 

and an additional holon-spinon {h/s) gauge invariance arising from the ambiguity in the 
decomposition ( p.6|) : 

h/s : Hje'^\ H* H*e-''^J, 

s,,-^s,,e*^^ s*,-^s*,e-<^ o e R. (2.12) 

It is important to remark that the theory in terms of {H, H*, S, S*, B, V} is equivalent 
to the original fermionic theory only if the h/s gauge is fixed to respect U{1) x SU{2) 
invariance. The h/s gauge-fixing will be discussed later. 

We first gauge-fix the U{1) symmetry imposing a Coulomb condition on B (from now 
on /i = 1, 2): 

d''B^ = 0. (2.13) 

To retain the bipartite lattice structure induced by the AF interactions, we gauge-fix 
the SU{2) symmetry by a "Neel gauge" condition: 

S,=crl^lQ, S* = (l,0)al^l (2.14) 

where |j| = ji + j2- Then we split the integration over V into an integration over a field 
V^^\ satisfying the Coulomb condition: 

Qt^yic) ^ 0, (2.15) 

and its gauge transformations expressed in terms of an S'f/(2)-valued scalar field g (not a 
second-quantized operator), i.e., Va = g^V^^^'g + g'^dag, a = 0, 1, 2. 
Integrating over Bq, we obtain 

B^ = B^ + 6B^, 6B^{x) = WH*H,d^aTg{x-j), (2.16) 

j 

— i f B 

where 5^ gives rise to a vr-flux phase, i.e., e ■''^p = — 1 for every plaquette p. 
Integrating over Vq, we find 

Vi'^ = E(l - H*H,){a\^\g]^g,a\^)ud,e.Tg (x - j)a^, (2.17) 
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where cTq, a = x,y,z are the Pauh matrices. After the x SU{2) field being gauge-fixed, 
as discussed above, the action ( p.lO| ) becomes 



S{H, H*, 



dx'{Y.[H*{do - S)h, + (1 - H*H,){al^\g]dog,al^^)n 

)9j(^x^)n + h.c. 



<ij> 



j<ij> 



(2.18) 



Here the boundary terms are omitted and the 5"^ x topology of the involved euclidean 
space-time imposes the vanishing of the topological term Tr JoxR^io^dg)^, where D is a disk 
of radius p. 

Equation ( |2.18| ) is the starting point for our subsequent analysis. In S{H, H* , g) the 
charge degrees of freedom (d.o.f.) are described by H,H* (2 d.o.f.) and the spin degrees 
of freedom by g (3 d.o.f.) subjected to a constraint (-1 d.o.f.) coming from the h/s gauge 
fixing, reproducing the correct counting of degrees of freedom of the original fermionic fields 
"^a, (2+2 d.o.f.) in the euclidean path-integral formalism. 



III. OPTIMIZATION OF THE SPINON CONFIGURATION 

To analyse ( p.l8| ) we first recall the strategy adopted for the ID case [^. We noticed 



that one can find an upper bound for the partition function of holons in a spinon back- 
ground. Moreover, one can find explicitly the spinon configuration, exactly saturating this 
bound. Then one can consider the quantum fiuctuations around this optimal configuration 
to evaluate related physical quantities. Here we will follow a similar strategy, namely to 
search first for the upper bound of the partition function for holons in a spinon background. 
It turns out that such an upper bound exists. However, unlike the ID case, we cannot find 
a spinon configuration exactly saturating this upper bound. Nevertheless, we can find a 
holon- dependent spinon configuration g"^ which is optimal, saturating the upper bound on 
average, and take it as the starting point to consider the spinon fiuctuations. 



A. Auxiliary Lattice Gauge Field and Upper Bound 



To find the optimal configuration we introduce an auxiliary lattice gauge field {A,U}, 
with Aj G R (real), f/<ij> G C (complex), |?7<ij>| < 1, and an action S = Si + S2, 

Si{H, H\A, U) 

= dx' {Y^[h; {do - 5) H,+ til -H*H,) A,] + J2i-tH:U^,,^H, + h.c.)}, 



S2{H,H\u) = j^dx' J2 ^{i-h*h;){i-h*h^){\u. 



<i3>\ 



(3.1) 



The action (|3.1| ) equals ( |2.18D if we make the identifications: 
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t/<.> - e-^/<^.>(^ + ^^)(al^l,t(Pe^/<.> ^'^)g,afU (3-2) 

(but in the derivation of the following bound these identifications are not made.) 
Let 

we prove in Appendix A the following upper bound: 

< j OTPiy*e~['^i^^'^*'0'^) + '^2(i^',i^*,0)]^ (3 3) 

where U is the time-independent [/-configuration maximizing 

/ P/JOT*e-[^i(^' 0' U) + S,{H^ H*, 0)] |^^^^^_ (3 4) 

To discuss the properties of U we first notice that the quantity optimized by U is the 
free energy F{U) of a gas of spinless holes at temperature T = f]~^ with chemical potential 
6 { 6 = corresponding to half-filling) on a lattice, with hopping parameter on the link 
< ij > given by t|?7<jj>|, in the presence of a constant (but not uniform) magnetic field 
where the flux through a plaquette p is given by arg (Ugp) and is subjected to an attractive 
n.n. density-density interaction with couphng constant j. We consider the system at large 
(3, and small 6 and make the following: 

Assumptions: 1) we consider negligible the density-density interaction, since the holon 
density (5) is small; 

2) we assume translational invariance of U. 



Remark: Assumption 2) appears to be reasonable in the light of the results of [^], where 
it has been shown that, the configuration U maximizing the determinant of the hopping 
matrix of the above system is translation invariant. 

By gauge-invariance (see Appendix A) the result of optimization depends only on |f/<ij>| 
and arg (Udp). As a consequence of the assumption 2), F{U) is monotonically increasing in 
|f/|, hence 

|f/<,,>| = 1. (3.5) 



It has been conjectured in |55] and proven in |56] that the ground state energy at T = 



of the system under consideration is optimized in the magnetic field chosen as: 

arg {Uqp) = 7r(l - 5) (3.6) 

which is the commensurability condition for the flux. It is then natural to conjecture that 
this remains true for large enough /3. 
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B. Optimal Spinon Configuration 



In the last subsection we have stated a kind of "theorem" (proven in Appendix A), now 
we will find out the consequences of this theorem in our context, namely, assuming ( p.5| ) and 
( p.6|) to hold, we attempt to find a holon-dependent spinon configuration {g"^) saturating 
the bound (pj.3D, using the identifications (|3.2|). 

Following a strategy developed in ID |^^, we introduce a first-quantized (Feynman- 
Kac) representation of the holon partition function in the presence of a (^-background. As 
in ID, a key ingredient in the analysis is that, for a fixed link < ij >, in every holon 

configuration the term {al^^glPe ■^<'J>^ (7jO"[?')ii appears, in the first quantized formalism, 
either in the worldlines of holons or in the Heisenberg term, but never simultaneously. This 
is the consequence of the single-occupancy constraint and permits a separate optimization 
of Si and 5*2, as required in the bound ( ^.3| ). It turns out, however, that, contrary to the 
ID results, we cannot identify a specific configuration g"^ saturating the bound exactly, but 
only approximatively, in an appropriate average sense. 

The whole procedure in 2D can then be justified in the limit t » J because the effective 
mass of holes is very heavy, as a result of the large number of soft spinon fiuctuations 
surrounding the hole in its motion |57] and in a sense our treatment can be considered as a 



kind of Born-Oppenheimer approximation for the spinous in the presence of the holons. 
Here we do not give details of the derivation, but rather introduce notations and quote 



the obtained results. Those interested in further discussions are referred to ||58| , |38| , |39| ,^ 
Let A denote the 2D lattice laplacian defined on a scalar lattice field / by 

j:\j-i\ = l 

let dii{uj) denote the measure on the random walks u) on the 2D lattice such that 



in. 



L(o)=i M^)^ > 0; 

let P/v be the group of permutations of elements and, for vr G P/v, let a^n) denote the 
number of exchanges in tt, then the partition functions of holons {H) in a given g background 
can be rewritten as: 

<ij>eui-^ 

where identifications ( ^.2| ) are understood and, for a fixed N, u = {ui, ...,un} denote the 
worldlines of holon particles, ur^ the components of a; perpendicular and u;" parallel to the 
time axis, respectively. 

To saturate the bound ( |3.3| ) with a configuration g'^{uj) we first impose 

iAj = {a\h]dogj(Tf)ii = Q, jeJ. (3.8) 
Eq. (|3.8|) is satisfied choosing g"^ constant during the period when no particle hops. 
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Imposing ?7<ij> = in 5*2 (see corresponds in the first quantized formalism to 

setting: 

(alSIPe*^<-> ""'^.al^l)!! = 0, < y > n a; = 0. (3.9) 

In physical terms this means that the s + id RVB order parameter [47| is very small (see the 
discussion at the end of Sec. IV). 
We notice that if 

gj = cos II + isinOja^, j ^ id. (3.10) 

for some angle 9j E [0,27r), ( PT^ ) is then satisfied. In fact, since V^'^^ depends only on sites 
where there are no holes ( see ( p.l7|) ), from (|3.10|) it follows that correspondingly 

V^'\x)=J2il-H;H,y—^ d.^Tg {x - j)a., (3.11) 
i ^ 

so that g\Pe^<^i>^ gj has only diagonal components. 

This result shows that in the Neel gauge, quite independently of doping concentration 
(since the condition f/ = in 5*2 does not depend on small doping assumption), one should 
expect that the physics is dominated by V^''^ only in the U{1) subgroup of SU{2) related to 
the axis chosen in the Neel gauge. 

The condition \tJ<ij>\ = 1 in 5*1 (see (|3.1|) ) corresponds to imposing: 

\{afg\Pe'^<^o> ^^"^ g.a^^U = 1, < ij >E (3.12) 



which means in physical terms that the AM order parameter is of the order 1 (see the 
discussion at the end of sec. IV). 

To discuss ( p.l2| ) we recall that the paths on which dfi^u) is defined are left-continuous 
so that at the jumping time r, co'j.(r) = lime\o uJriT + ^), or, in simpler terms, one should 



think the holon at r at the end of the jumping link, oriented according to the increasing 
worldline time of the holon. As a consequence, if < ij >E ur^, either i E u or j E u, but 
never both. Let us assume j E then according to the previous requirements 

gjPe ■^<*i> ^^'^ = cos O^ij^l + i sin 9<ij>a„ 

for some angle ^<ij> E [0,27r). We represent 

gj = cos ipj + ia ■ rij sin Lpj, 

for some unit vector fij and angle ^pj E [0, 27r). 
From ( p.l2| ) one immediately obtains 

vr 

V'i = 2' = 0. (3.13) 



Finally, let us try to impose the condition ( |3.6| ). This translation invariant condition 
cannot be exactly fulfilled (unlike in ID where a configuration g^ exactly saturating the 
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bound analogous to (|3.3|) can be found). However, we notice that the -B-dependent part of 
arg (Ugp) has a translation invariant mean satisfying the above condition, ( see (|3.3| ) and 
( p.l6| )), hence it is natural to impose (on average): 



(c7l'l^le'^<«J>^^'^(7,al^l)^^~l, <zj>eui-^. (3.14) 

Defining 

-f5:(-l)V,arg {i-j) 
9j = e , (3.15) 



and choosing 



^, = r^' (3.16) 



we can kill the fast fluctuating first term in ( p.ll|) . The remaining term, denoted by 



'-1)1^-1 



j ^ 

for small hole concentration is a slowly varying field yielding a contribution 0[S) to 
arg(?7<ij>), with zero translational average. The final result is that we can assume for 
the optimizing configuration, using ( |3.13|) , ( |3.15| ), ( |3.16| ) 



and we immediately derive from (|3.14| ) the condition Uj^ = 0,njy = (— l)!-'!. 



Notice that from these definitions: 

^,4^1 = (3.17) 

In view of the optimization discussed above it appears natural to introduce a variable 
Rj G SU (2) describing spinon fluctuations around the optimizing configuration, through the 
definition 

-|y(-l)V,arg(£-j) 
g, = -g,Rfg, = e i?/ 2 (-1)'^'^.^. (3.I8) 

with Rj being represented in CP^ form as 

where bj^ is a two-component complex field. The optimal configuration g"^ is given by 
R = l. 

Using eq. ( p.l7| ) and the S'[/(2)-gauge invariance of the (formal) measure Vg to absorb 



g, the partition function of the t — J model can be exactly rewritten in terms of the euclidean 
action S = Sh + Sg, 
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+ ^ [^tH;e-'^<^J>^^ + ^^^H,{a]i\RlPe'^<v>(^ + ^^^R,cT\^\)n + h.c.]}, (3.20) 



<'j> 



+ E ^{1-H:H,){1-H*H,) [|(al^li?JPe^^<^^>(^ + ^^)i?,ai?l)np-y}, (3.21) 



<*j> ^ 



where 6V = V^"^ - V. 

Let us remark again that in ( p.2(]| ), (p.21|) no approximations have been made. The point 



of the above analysis is that for large enough (3 and small enough 6 we expect that R and 
6V describe small fluctuations. 



IV. SPINON EFFECTIVE ACTION 
A. The Main Approximation 

To proceed further we make the following 

Appr oximat ion {6V = 6B = 0) : we assume that the spinon fluctuations (i?) are small 
enough (for P large, 6, J / 1 small) that we can neglect their back reaction on the gauge field 
V, i.e. we set 6V = 0. We expect that the main effect of the neglected fluctuations of V is 
to convert the gauge invariant spinon field operator reconstructed from the euclidean field 

into a semion field operator. Retaining the fermionic nature of is then inconsistent with 
neglecting 6V (which also introduces a fictitious parity breaking) unless we neglect also SB, 
responsible for the semionic nature of the gauge invariant holon field operator reconstructed 
from the euclidean field 

In ID the proper account of the statistics of the holon and spinon field operators was 
crucial for deriving ( within the C-S approach) the correct physical properties of the model, 
known by Luttinger liquid and conformal field theory techniques |^5|. However, in 2D we 



believe the statistics of ho Ions and spinous is less crucial because we expect that, contrary 
to ID, they form a bound state, as will be discussed later. 

To derive the low-energy spinon action let us start computing the link variable 
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where a<ij> = e2 J<«i> ^ . Looking back at equation ( p.21| ) we find that in the hopping term 



of holons only the diagonal elements of (|4.1|) appear, a kind of gauge-invariant AM variable 
J^, whereas in the Heisenberg term only the off-diagonal elements of ( [4.1D appear, a kind 
of gauge-invariant RVB variable 1^^. According to the optimization arguments given in the 
previous Section, the vacuum expectation value of the AM gauge variable is expected to be 



s-like, real and close to 1 ( see eq. (|3.14| )), while the RVB order parameter should be rather 



small( see eq. (|3.9| )). These anticipations are fully confirmed by the mean field calculations 



B. Nonlinear a Model with Mass Term 

We now derive a low-energy continuum effective action for spinous by rescaling the model 
to a lattice spacing e << 1 and neglecting higher order terms in e. As it is standard in AF 
systems [^, we assume 

h)J^ph,pr^Vl, + {-l)\^\eL,, (4.2) 

with fi^ = / < 1, ■ L = 0, where f2,L are defined on a sublattice e.g. Vtj = i^j^+i j^) 
Lj = Lj_^_^_i j^, ji = j2 mod 2 and they describe the AF and ferromagnetic fluctuations, 
respectively. It is useful to rewrite f2 in the CP^ form: 

— * 

^ = Z*aal3Zf3, Z*^Za = /, (4.3) 

with a = 1, 2 as a spin ^ complex (hard-care) boson field. Consistently with the slowly- 
varying nature of V for small hole concentration, we assume 

e-'I<ji>y^ _ 1 ^ e(-zK)(j) + '-{-^V,f{J) + 0{e'). (4.4) 
On the rescaled lattice the Heisenberg term becomes 

= ^ E + 2L^ + V^U) [(fi,.)^ + (fi,,)^] I + 0(e). (4.5) 

For the temporal term we obtain analogously 

- {a\^R]d,R,af)^, = Y.{-iy'^z*,doz,^ + ■ {Q, A doQ,) + 0{e'). (4.6) 

j ^ 

The first term in (E]B) actually vanishes, since otherwise it would produce a topological 



6'-term which is known to be absent in 2D 

In treating the holon density terms in the spinon action in MFA, we keep only the leading 
fluctuation terms and neglect terms of the order of 0(J5),0(5^). Then, integrating out L 
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and taking the continuum limit, from ( |4.2| ) - ( [4.6| ) we obtain the NLcx-model effective action 
for spinous S* + S'^, 



s: = f d'x-iidony + v',id,n)' + (nmn s: = -- f d'xnlv,', (4.7) 

•'9 9 



where the couphug constant g and the spin wave velocity Vg are easily derived functions of 
J, 6, e. 

We now make the Approximation P: we treat the term S'^ as a perturbation. To under- 
stand the physics described by S* we first notice that if were absent, the NLa-model 
would be in the symmetry broken phase, since g is small (~ J) at the lattice scale. For larger 
scales geff then flows towards its critical value, describing the large-distance properties of 
a NLcr model with an insulating "Neel" ground state and spin-wave Goldstone excitations 
To get an idea of the effect of we replace the NLcr model constraint = / by 



a softened version, adding to the lagrangian a term X{p? — fY, expected to produce the 
same low-energy behaviour, and we replace V^, a function of the holon positions, by its 
statistical average < >. By its definition, is positive definite and we give a rough 
estimate of < > by first performing a translational average at a fixed time over a fixed 
holon configuration and then an average over holon configurations with mean holon density 
5. 

In the first step let x = {xi}^i denote the holon positions and let q = {qi}iLi,qi = 
g^j^^ restrict the computation to a finite volme |A| and lattice cutoff e. Let 
arg*^, 9^, A*^ denote the angle-function, the derivative in the /i— direction and the laplacian 
in the e— lattice, respectively, then using the equality 

e,,dtarg'ix-y) = dliAT\x-y), 

we immediately obtain 

/ d'xV^\x) = E qMAr\x,,Xk) ^ E - x,\ + e). (4.8) 

|A| 7a |A| |A| 27r 

XCA £CA 

Equation ( [4.8| ) looks like the energy per unit volume En^a{x) of a neutral two-component 
system of particles with charges ^ volume |A| interacting via 2D Coulomb 

potential with ultraviolet cutoff e. The average of EN,A(yX) over the positions of the charges 
in the limit |A| y R^, N y oo with fixed average density 6 can be identified with the free 
energy of the above system at /5 = 1 and chemical potential S in the thermodynamic limit, 
and this gives our estimate of < >. The behaviour of the free energy can be obtained 
through a sine-Gordon transformation 0] and for small 6 it is given by ||61|| : 

<V^>^-6\n6. (4.9) 

Hence in our crude approximation < > acts as a mass term increasing with 6. If we 
assume that the scaling limit and perturbation expansion in < > commute with each 
other, we find that its effect is to drive the NLcr model at large scale to the disordered massive 
regime with a mass gap for Q of the order of rri^iS) ~ —5 In 5. Within this approximation, 
for a more careful analysis one should consider the renormalization group equations for g 
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and A including the perturbation < >-term from the beginning. Then the value of m^((5) 
could be renormalized. 

A slightly better approximation is to consider the holons as slowly moving randomly 
distributed impurities, creating a random potential and analyse the behaviour of the spin- 
wave Q in the presence of this potential (this can again be justified for the limit J/t << 1 
with a large effective mass of holes). At a fixed time the random potential behaves as 



\x — \y — w\ 



Hence it is positive and roughly falling like r ^, where r is the distance to the closest 
impurity. This kind of systems have been considered in and it is expected to be in the 



localized regime for f2, where (random) averaged Green functions exhibit a mass gap m{6) 
roughy proportional to the inverse mean free path. Hence up to logarithms m1{5) ~ 5 for 
small 8. 



C. Crossover from Short Range to Long Range Antiferromagnetic Order 

All the above arguments strongly suggest that the spinon system described by the action 
S** in ( [4.7|) exhibits a mass gap ms{5) increasing with 5 (at least for sufficiently large, 
but, nevertheless 5 « 1), vanishing at 5 = 0, thus showing an expected crossover to the 
insulating "Neel state" at zero doping As a result of m((5) > 0, the constraint VL^ = f is 



relaxed at large scales ||2^ and a lagrange multiplier field introduced in the action to impose 
the constraint would just mediate short-range attractive interaction between spinous with 
opposite spin, with strength 0(ms(5)~^) rotationally symmetric in spin space (see for 
the discussion of an analogous situation). Another short-range interaction, but uniaxial in 
the spin space, is introduced by the perturbation term S'^ in (|4.7|) , with strength Q{ml{5) / g) . 

We can summarize the above discussion by rewriting the NLa model effective action for 
spinous in the CP^ form, neglecting short range interactions, as 



J[0,/3]xR2 g 



|(9o - zld,zp)z^\' + vl\{d^, - zld^zp)z^\^ + ml{6)zlz^ . (4.10) 



In the NLcr model without mass term {6 = 0) the constraint z^Za = f and the symmetry 
breaking condition, e.g. < Z2 >7^ 0, lead to excitations described by a complex massless 
field S" ~ < ^2 > with massless "relativistic" dispersion relations, corresponding to spin 
waves. In the NLa model with mass term the absence of symmetry breaking and the effective 
softening of the constraint lead to excitations described by the spin ^ two-component field 
Za, with massive "relativistic" dispersion relations. However, the self-generated gauge field 
z^daZa,a = 0, 1,2 would confine the spin | degrees of freedom into composite spin 1 spin- 
wave fields. On the other hand, as we shall see later, the coupling to the holons will induce 
deconfinement of spin ^ excitations. 

To give a kind of "microscopic" interpretation of the above results using the slave fermion 
picture in terms of the hard-core bosonic fields ha,h*^, we follow the approach of Yoshioka 
p!3| . First consider the case 5 = 0; then a MF treatment with an s-like RVB order parameter 



yields an energy gap vanishing at (|, |), (— |, — |) in the Brillouin zone and an expansion for 
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low momenta around these 2 points shows that the corresponding excitations are described 
by a complex massless field with relativistic dispersion relations, which can be identified with 
a spin-wave field S; the ground state is then the insulating Neel state . A MF treatment 
with an s + id RVB order parameter yields an energy gap vanishing at (±|,±|) in the 
Brillouin zone, and an expansion for low momenta around these 4 points shows that the 
corresponding excitations are described by a two-component massless field with relativistic 
dispersion relation which can be identified with a spinon field Za [0. These excitations 
are turned into massive ones in our approach by the ms{6) term. In some sense, our result 
in the NLo" model corresponds to a kind of slave-fermion approach with a gauge-invariant 
RVB order parameter s + id-like and minimal gap ^^(5) at (±|, ±f )• We should emphasize, 
however, that our approach is different from the slave-fermion formalism which considers 
long-range ordered AF state with gapless excitations, while we consider short-range ordered 
AF state with gapful excitations. 

To the best of our knowledge, the present formulation is the first successful attempt to 
explicitly include the AF fiuctuations into the RVB-type scheme of treating the strong corre- 
lation effects in a self-consistent manner. Upon doping the long range AF is destroyed, being 
replaced by short-ranged AF order, which is physically obvious. We have obtained an ex- 
plicit doping dependence of the gap value which has the correct extrapolation at zero doping 
(gap vanishes) [Q. The doping dependence of the AF correlation length ^ ~ ~ 
expected from our calculation, agrees very well with the neutron scattering data Here 
we also propose a new interpretation of the spin gap effect in underdoped superconductors 
- mainly due to the short range AF order. The recent numerical simulations ||68| seem to 
support our interpretation. 



V. HOLON EFFECTIVE ACTION 
A. Holon-Spinon Coupling via Self-Generated Gauge Field 

— i f B . 

Now turn to holons. We use a gauge choice of the field B such that e j<»j> is purely 
imaginary and assume, following eq. ( p.l4| ) - ( p.l7|) that the gauge - invariant AM parameter 



is ~ 1, which permits us to use the equality (see eqs. (|4.7|), ( p.l9| )) 

{a<ij>biibji + «<ij>&»2&i2) ^ 1 = b*^bii + 6*2&i2. (5.1) 

We obtain a low energy continuum effective action for holons by rescaling the lattice 
spacing to e << 1 and neglecting higher order terms in e. Making use of ( p.l|) one obtains 



<'j> 

'bja bia , /b^^ ^j<^\l, ]^^^^\ 



+ {H:H, + H*H,) [(&L(^^^) - I + 0(e) I , (5.2) 

where #(j) denotes complex conjugation if j is even. Neglecting the 6-terms, the action 
( ^.2| ) describes the usual two-component Dirac ("staggered") fermions of the flux phase, 
with vertices of the double-cone dispersion relations in the reduced Brillouin zone, centered 
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(in the e = 1 lattice) at (±f , ±f ) and chemical potential S. To define a continuum effective 
action, using the standard procedure for the fiux phase |2^, we first define 4 sublattices: 1) 
for ji, j2 even, (2) for ji odd, j2 even, (3) for ji even j2 odd, (4) for ji, j2 odd. They can 
be grouped into two "Neel" sublattices A = {(1), (4)}, 5 = {(2), (3)}. The holon field H 
restricted to the sublattice # is denoted by H^^\ The holon action can then be written in 
the continuum limit as a bilinear form in H* and H with kernel: 



f do- z^doZa - 5, 
it{di - z^diZa), 
-it{d2 - zld2Za), 
\ 0, 



it{di + z^diZa), -it{d2 + zld2Za) 

do + z^doZa -S, 0. 

0, do + z^doZa - S, 

lt{d2 + Z*^d2Za,), lt{di+Z*^diZo^), 



\ 

lt{d2 - Z^d2Za) 

it{di - z^diZo), 
do - z^doZa - 6 j 



(5.3) 



where we used the decomposition of ha carried out in the previous Section (all ferromagnetic 
terms are 0(e) and so they do not appear in the continuum limit (|5.3|) ). 
Setting 



7o 



^1 



7/. 



(A) 



(A) 



[ay, a. 



4 = 7m^/., 

e-fi7(3) + e*fi7(2) ) 

e-fi/(2) + e*f/7(3)x 
e-fi7(4) + e*fi7(i) ) 



# 



7o*# 



(5.4) 



and assigning charge ca = +1,6^ = —1 to the fields corresponding to A and B sublattices, 
respectively, the continuum effective action for holons can be rewritten as ||69|| : 



St 



[0,/3]xR2 



7o(5o 



6 - CrZl^doZ, 



QZq^ 



+ +t{^ - erZ*J}Za) 



(5.5) 



Hence, describes the coupling of the holon field to the spinon-generated gauge field 
z*JiaZa in a "relativistic" Dirac - like form, with opposite coupling for the two Neel sub- 
lattices. From ( |5.5| ) it is clear that (neglecting the gauge terms), the upper components of 
describe gapless excitations with small FS [ep ~ 0(5t)), whereas the lower components 
describe massive excitations with a doping dependent gap mh{S) ~ 2\/S. In terms of lattice 
fields we have excitations of the two bands (gapless and gapful) supported in the reduced 
Brillouin zone due to the presence of the Neel sublattices. This leads to a band-mixing due 
to the non-diagonal structure of the 7-matrices and this in turn yields a "shadow band" 
effect, i.e. a reduction of the spectral weight for the outer part of the F.S. facing (vr, vr), 
in agreement with experimental data |]10[ and MF numerical simulation [Q. The FS of 



underdoped superconductors has also been considered in the SU (2) gauge field theory p2| . 
However, there are some differences between the results obtained there and the predictions 
of the present model (for details see @). 
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B. U{1) X SU{2) Gauge Invariance of the Total Action 



Neglecting the quartic fermion term, which would produce a short-range repulsion be- 
tween holons with opposite charge, the full continuum effective action can be rewritten 
introducing an auxiliary U{1) gauge field A (over which one integrates in the path-integral): 

+ml{6)z*^z^] + E'=i ^rboido - erAo -6)+ - e,^)]^,.} (5.6) 
The action S* is invariant under the U{1)- gauge transformation 

v|/,(x)^e^^-^(^)vl/(x), 
^,.(a;)^e-^^-^(^)^(x), 

Aa{x) Aa{x) + daA{x), A(x) G R. (5.7) 

Going backwards we recognize in ( ^.7|) the continuum limit of the h/s gauge invariance 
(^TT^) which is gauge-fixed imposing, e.g., a Coulomb gauge condition on Aa- 



Remark - One easily finds that the lattice counterpart of the A^ gauge-fixing is a gauge- 
fixing for arg (S*Pe*'^<*^> ^Sj) which is U{1) x S'f/(2)-gauge invariant as required in Remark 
in Section II. 

In fact 

arg(E*Pe'-^<*J> ^S,) - arg {a<ij>b*,bji + a*^,jyh*^hj2) ^ z^d^z^ - A^, 

where in the last equality we have neglected a quadratic fermion term fictitiously introduced 
in S* (compare it with ( |5.5| )). 

We can use S* to compute, as in EO], the gauge field propagator induced by the spinon 



and holon vacuum polarizations. In the Coulomb gauge, denoting by 11 (n" ) the transverse 
(longitudinal) polarizations, respectively, we have 

{A,{u;,q)A,{-u;,-q))-'^ ^^^^^ ^ 



{Ao{uj,q)Ao{-uj,-q)) 



g2 ; {Uj + Ui){u,qy 



q^{ul + ul){u,qy 

Since the spinon is massive the vacuum polarization of the spinon system alone would 
be Maxwell-like in the absence of holons and it would induce a logarithmic confinement of 
spinous (n-*- ~ o;^ + g^, n" ~ cu^ for small to, q). However, due to their gapless component 



the holons induce a transverse vacuum polarization with Reizer singularity |21| (11^ 



1^1 + X^f') ci-nd a leading behavior in Ilj^ ~ ^ for small uj/q, leading to a short range Aq 
propagator, so that the full gauge interaction is not confining |7^. 

As a result, holons and spinous are true dynamical degrees of freedom in the model. 
Nevertheless, since we are in 2D, the attractive force mediated by the gauge field is expected 
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to produce bound states with the quantum numbers of the spin wave (see for the 
discussion of a similar problem) and presumably (see e.g. [^) of the electron. This could 
explain why neglecting the (expected) semionic nature of holons and spinous {6V = = 6B) 
can be justified to some extent in 2D. On the contrary, in ID the spinous and holons are 
deconfined, and their statistical properties as semions are crucial for a proper account of the 
spin-charge separation. 



VI. CONCLUSIONS 

To summarize our results it is interesting to compare the features appearing in the 
(underdoped) 2D and ID t — J model (for small 6,J/t) within the U{1) x SU{2) gauge 



approach followed here and in K 



In both cases, applying to the t — J model the U{1) x SU{2) C-S bosonization in terms 
of gauge fields B,V, we separate spin and charge degrees of freedom of the electron, and 
describe them in terms of spinon (z) and holon (\1/) fields. 

1) The low-energy properties of the spinous can be described by a NLcr-model: in 2D it 
is massive, and the mass gap (vanishing at 5 = 0) is due to the coupling to holons mediated 
by a non-pure gauge SU (2) field V; in ID no such non-trivial gauge field exists and the 
spinous are massless. The self-generated spinon-gauge field {z^daZ^) would confine the 
spinous into spin waves, but this is prevented in 2D by the coupling to holons, while in ID 
by the presence of a topological 6 = n term in the NLa model, which is absent in 2D. 

2) The low-energy properties of holons are described in 2D by a Dirac action, inducing 
band mixing, whose appearance is due to the presence of a non-pure gauge field B (charac- 
teristic of the flux phase); the absence of such a field in ID implies that holon are described 
in term of a spinless fermion action. 

3) Spinous and holons in 2D are coupled by the spinon-generated gauge field, carrying one 
degree of freedom (after gauge fixing) and this presumably induces binding; the gauge fixing 
kills the degree of freedom of an analoguous gauge field in ID and as a conseguence, spinous 
and holons are free. The statistics is then dictated by the U{1) x SU{2) decomposition and 
the corresponding (gauge-invariant) fields obey semionic statistics. 

To reiterate, we outline some of the distinct features of our present study: 

- The short range AF order is the main reason of the existence of the spinon gap (at 
least in the underdoped samples) which can explain a variety of experimental observations. 

- Neglecting gauge fluctuations the 2-point correlation function for the electron exhibits 
half-pocket FS in the reduced Brillouin zone, with a pseudogap minimal on the diagonal, 
"shadow band" effect and a quasi-particle peak due to the holons. These features have been 
demonstrated in the MF computation of , where a similar action for the system has been 



used, with a twist in holon-spinon statistics, i.e., the spinous are fermions, while the holons 
are bosons. 

- The spin-wave persists even in the region without AF long range order, but as a short 
ranged field. 

- We expect (this is still under investigation) that the "composite electron" exhibits a 
Non-Fermi liquid behaviour due to the appearance of Reizer singularity, thanks to transverse 
gauge fluctuations, in the self-energy of its holon constituent. 
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To conclude, the C-S bosonization approach, in spite of its technical complications, is 
promising in providing a natural framework to describe the normal state properties of under- 
doped superconductors. After having settled in this paper our framework for a discussion 
of the t — J model at small 5 and J/t as a model for underdoped high cuprates, in 
the forthcoming paper we will compute correlation functions and compare the results with 



experimental data [50 
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APPENDIX A: 

The proof of the bound (|3.3D uses techniques adapted from the proof of the diamagnetic 



inequality given in [72 



The main ingredients are: 

1) invariance of the action (|3.1| ) under the gauge transformations: 

A, ^ + ^oA,-, ^ t/<.,>e^^^* - , 

Hj^Hje^^i, H*^H*e~^^i, A^ G R. (Al) 

2) reflection (O.S.) positivity in the absence of gauge fields: define an antilinear involution 
9 on the holon fields supported in the positive-time lattice (the time interval is identified 
with [-f,f])by 

9:H,{r)^H';{-r), H^^r) ^ H,{-r) 

d{AB) = 9{B)9{A), A,B ej+ (A2) 

where J7+ denotes the set of functions of H, H* with support on positive time. Reflection 
positivity is the following statement: VF G J7+ 

{F9F) = J VHVH*e'^^^'"*^F9F* > 0, (A3) 

where S{H, H*) = S{H, H*, 0, 1) (see eq. (PD). 

The involution 9 defines an inner product (■, ■) on J'^ by {■9-) and as a consequence, 
Schwartz inequality holds: WF, G E J+ 

\ {F9G) \ < {F9F) 5 {G9G) ^ . (A4) 



Although the proof of (^) is standard [^ , let us sketch it for readers' convenience and 
to set up a notation useful in the proof of ( p.3| ). 
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We discretize the time interval introducing a time lattice with spacing e = f32~^ and we 
choose the time— plane lying halfway between lattice planes. We denote by S^^S'L and 
S'^_ - the terms of the discretized action existing in the positive, negative time lattice and 
coupling the two, respectively, then we have: 

^1 = esxs%^ = -Y.[H-{{)H,{^) + H,{^)H*{^] 
j 

= - E [H;C^)eH;C-) + H,C-)9H,C-)\ . (A5) 

j 

Let P+ denote the set of functions of H, H* given by linear combinations with positive 
coefficients of elements of the form F6F, FtJ^. It is easy to see that P+ is closed under 
multiplication and summation with positive coefficients and for CiFiOFi G P+ we have 

j {VHVH*),Y.('iP^^F^ = ^''^l / ('^HVH*),F,\^ > 0, 

i i 

where [VHVH*)^ denotes the (formal) measure on H, H* in the path-integral for the dis- 
cretized time model. 

From (^) it follows that, for F e J+, e'^"^"^"*^ FOF e P+, so that 



J {VHVH*),e-^'^"'"'^FeF > 



and taking the limit e \ we obtain 

We turn now to the proof of ( ^.31) . We work again with the discretized time lattice. 
First we use the gauge invariance (|Al|) to set A = on the links in 5"+-. We denote by 
{A^^\u[^^} {{A^\U2'^^}) the restriction of the gauge fields to the positive (negative) time 
lattice and we set 

^(1) = ^-[SUH,H*,a['\u[''>) - SUH,H*)] 

where r is reflection w.r.t. the time-0 plane. Applying ( [A4| ) we derive the upper bound: 

\E'{A,U)\ = I J OTOT*e~'^'(^'^*'^'^)| = 

= \{F^'^eFi'^)\ < (P«eP«)i(P«^P«)i 

We now choose a time plane, different from the time— plane, lying halfway between 
lattice planes, use the time translational invariance of S (inherited from the antiperiodicity 
in time of H, H*) to bring this plane to the position of the time— plane shifting the whole 
lattice and repeat the above procedure both for Fi^'^OFi^^ and F2^^9F2 (which are no more 
in the form F9F with 9 referred to the new time— plane). As a result, we obtain, with 
obvious notations, Ff \ i = 1, 2^, and then iterate. Finally we derive the bound 
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2N 

mA,u))\<i[{FroFrr''- (a?) 

A close inspection shows that the gauge fields appearing in F^^^OF^^^ have the following 
properties: A = and U is time independent. Using a hamiltonian transcription, with 
obvious meaning of notations, looking at we improve the bound by means of Golden- 

Thompson inequality: 

(iT^W^iT^W) = Tre-/^{[^i + = 0)] + [H^ - H^iU = 0)]}^-(3[H2 - H^iJJ = 0)] 

< T^^-mi + H2{U = Q)]^-m2 - H2{U = 0)] < Tre-^'yH^ + ^^(f/ = 0)]^ (A8) 

where we used [H2 — H2{U = 0)] > 0. Eq. ( [A 81 ) is equivalent in the path-integral formalism 
to bounding the r.h.s. of ( |A7| ) setting everywhere U = in 82- Let us define 

E'{U) = J {VHVH*),e-^i(^^ 0' ^) + ^'^(^^ 0) la^t/^o (A9) 

and let U denote the gauge field configuration maximizing 'E^{U), then the r.h.s. of ( |A7D , 
using ([A8|) , is bounded by S^(f/) and taking the limit e \ we recover (p.3|). 
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